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Abstract 

The one loop corrected effective Lagrangian for the quark-squark-chargino and quark- 
squark-neutralino interactions is computed. The effective Lagrangian takes into account 
the loop corrections arising from the exchange of the gluinos, charginos, neutralinos, W, 
Z, the charged Higgs, and the neutral Higgs. We further analyze the squark decays into 
charginos and neutralinos and discuss the effect of the loop corrections on them. The 
analysis takes into account CP phases in the soft parameters. It is found that the loop 
corrections to the stop decay widths into chargino and neutralinos can be as much as 
thirty percent or even larger. Further, the stop decay widths show a strong dependence 
on the CP phases. These results are of relevance in the precision predictions of squark 
decays in the context of specific models of soft breaking in supergravity and string based 
models. 



^: Permanent address of T.I. 



1 Introduction 



The qq'iXj^ and qqix'j interactions are of great interest since they enter in the decay of 
squarks. We expect that such decays will be observed at the coUider experiments. Specif- 
ically one expects under the usual naturalness criteria that most of the sparticles should 
become visible at the Large Hadron Collider (LHC) with the possibility that some of the 
sparticles may also become visible at RUN 11 of the Tevatron. Measurements of sparti- 
cle masses and of their decay branching ratios will be a primary focus of attention after 
the discovery of such particles, while a more precise measurement will come eventually 
at the next linear collider (NLC). With the above in mind it is of great importance to 
refine the theoretical computations of the decay branching ratios beyond the tree level 
predictions. In this paper we extend the previous analyzes of squark decays into charginos 
and neutralinos[l, 2, 3] by taking into account the loop corrections with CP phases. For 
this purpose it is found advantageous to compute the one loop corrected effective La- 
grangian for the qq^Xj^ Q^Qi^j couplings. In the analysis we also include the effect of 
CP phases. It is now well known that large CP phases can be made compatible[4, 5, 6, 7] 
with the experimental constraints on the electric dipole moments of the electron[8], of the 
neutron[9], and of the Hg^^^ atom[10]. Further, if the phases are large they could affect 
a whole host of low energy phenomena. These include the effect on the higgs masses, 
couplings and decays[ll, 12, 13, 14, 15, 16, 17, 18], dark matter [19, 20], and a variety 
of other phenomena[21]. The outhne of the rest of the paper is as follows: In Sec. 2 we 
compute the effective Lagrangian for the thiXj and hiiX^ interactions. In Sec. 3 a similar 
analysis is done for the qqiX^j interaction. In Sec.4 we give an analysis of the decay widths 
of the squarks into charginos and neutralinos using the effective Lagrangian. In Sec.5 we 
give a numerical analysis of the size of the loop effects on the decay widths. We also 
study in this section the effect of CP phases on the decay widths. Conclusions are given 
in Sec. 6. 
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2 Effective Lagrangian for qq'jxf Interaction 

In this section we study the effect of loop corrections on thiXj and on hiix'^j interactions. 
We begin with the tree level Lagrangian density 

C = gt{RbijPR + LujPL)Xj^i 
+gb{RujPR + Lu,PL)Xjii + H.c. (1) 



where 



Rbij — —{UjiDhii — Ki,Uj2Dh2i) 

Lbij — KtV*2Dbii 
Rtij — —{VjiDtii — KtVj2Dt2i) 

Luj = KbU;,Dm (2) 



and where 



Ktib) = -^-^^ (3) 

y2mw sin (3 {cos (3) 

and the matrices f/, V and -D6(t) are the diagonalizing matrices of the chargino and squark 
mass matrices so that 

U*M^+V^^ = diag{m^+,m^+) 

D\MlD, = diag{ml,ml) (4) 

where m + (i=l,2) are the eigen values of the chargino mass matrix and m| (i=l,2) are 
the eigen values of the squark mass^ matrix. The loop corrections produce shift in the 
couplings of Eq. (2) as follows 

Ceff = gt{{Ruj + ARbij)PR + {Uij + AUij)PL)xpi 

+gb{{Rtij + ARtij)PR + {Ltij + ALtij)PL)xfi + H.c. (5) 

where AR^ij, ALinj, ARtij, and ALuj are the corrections that arise from the diagrams in 
Figs. (1-4). As is conventional we will use the zero external momentum approximation in 
the analysis of these corrections (see, e.g., Ref.[22]). 
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2.1 ARbij and AL^ij analysis 

Contributions to ARi^j and AL^jj arise from the nine loop diagrams of Fig.(l). We discuss 
now in detail the contribution of each of these diagrams, figs, (la-li). We begin with the 
loop diagram of Fig. (la) which contributes the following to ARbij and AL^ij 



k=l 



"^^ k=i 



where 



1 z X y 

f{x, y, z) = Azxln- + xyln- + yzln-) (8) 

[x - y){x - z)[z - y) x y z 



Next for the loop Fig. (lb) we find 



fe=i 1=1 



(/3«Aife + alDt2k)-^nml mjo, m|) (9) 

k=l 1=1 

rribm^o 

{atiDtik - ltiDt2k) ^g^a' /(r»b, m^o,m^J (10) 

^^^^'^ 2mw cos /9 (sin /9) 
/36(t)fe = e(56(t)-^ife H -^2U^36(t) - Q6(t) sin^ 6w) 

cos f7VF 

lh{t)k - eQb(t)Xik — Z ^2k U^J 

cos ny[r 

where X"s are given by 

= Xifc cos 9w + sin 
= -^ife sin 6*1^ + X2k cos 6*^^ (12) 
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and where X is the matrix that diagonahzes the neutrahno mass matrix so that 

X'^M^oX — diag{m^o,m^o, m^o, m^o) (13) 
Fig.(lc) contributes the following 

k=l 1=1 

iGki{Yi2 + iYis cos /3) + GUY12 - iYis cos /5) + Hf,i{Yn + iYjs sin /3) 

+H:,{Yn - iY,, sin /3) ) (C| + iC^) ^/(m?, mf^ , m% ) (14) 

where F is the diagonalizing matrix of the Higgs mass^ matrix 

YMl,ggsY^ = diag(ml^,mjj^,mjjj (15) 



and 



^ J^^n ((~^ + 1 ^w^)-D6ii-D6ij - J sin^ 9wDl2iDb2j) sin/3 
v2cos&^ 2d d 

+ ^ /^i^:H^.2, (16) 

V 2"7.vy COS p 

^ — J^^o + ^ ^w)DluDbij - \ sin^ ^^yL)*2i^fe2i] cos/? 

; [-'-'feli-^blj + ^b2i-^b2j\ 7= o^b2i^blj (-1- ' j 



V 2m COS /3 V im^y cos /? 



and 



Ctt - cos xt - C^i sin 
Ctf = sin + C^i cos Xt 
^2^1 = Re{ht + 5/it)Yi2 + [-/m(/it + 5ht) cos (3 
+Im{Aht) sin P]Yis + Re{Aht)Yn 
V2C^ = -Im(ht + 5ht)Yi2 + [-Re(ht + 5ht) cos /3 

+Re{Aht) smP]Yi, - Im{Aht)Yn (18) 



with 



7m(f^ + f^cot/3) 
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and 



h = '"^^ (20) 
\/2mw sin (3 



The corrections /Sht and 5ht are defined in Appendix A. 

ALg] = E E ii'tV^/2Aiik[G'fe^(r,2 + lYi^ COS /?) + GUYi-, - lY,-, cos /?) 

V2 

+zF,3sin/5) + i74(Fa -^Il3sin/?)](C| -^C,^)^/(m^m?^,m^^ (21) 
Fig. (Id) gives the following contributions 

= — E(/?6/Ah + ocliDm^{Bl: - Bl;)^,^ sin/3^^^/(m^, m^o, m^,.) (22) 



where 



Bl^ = -^(/ift + Shbje-'^'* sin /? + ^A/ifeC"^'^" cos /? 

-^(/it + 57^)6*^" COS ^ + ^Mte*^'* sin /3 
= + ^)e^^" cos /3 + ^M|e'^" sin /3 

+ ^(/i6 + ^)e-^^" sin/3 - ^A^e"^^"* cos/3 (23) 

where = (xb + Xt)/2 and where Xb is defined by the following 

7m(f^ + 4^tan/3) 

tan Xb = \t ^ (24) 

l + i?e(|f + f^ tan/3) ^ ^ 

and 

where the corrections A.hf, Shf, Ahf and 5hf are defined in Appendix A. 



where 



a/2 JL CD rribm^o „ „ „ 

= y E(«wAh - 76/^20(^6? + B^nCij cosP-^f{ml mjo, m|-) (26) 
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Next we discuss the contributions from Fig.(le). Here on using the properties of the 
projection operators, i.e., ■j^Pr = PlJ^, PlPr = 0, and the property of Dirac 7^* that 
QtivYY = 4, we get 

Ai?g = (28) 

and 

ALg = -Ag p R'ljiauDm - luD,2iy^f{ml mjo, m^vK-) (29) 

where 

R'ij — -^j^ziVjz + X-i^^Y (30) 

Contributions from Fig. (If) are as follows 

2 3 

= E Ef^iK^i + ilfes sin /5) + S^i{Yki + iYhz cos /3)] 

l=\ k=l 

{Cg + iCg)[UaD,u - K,UiAi]'^^^^^f{ml m\-,m\) (31) 



and 



ALg = ^ E E [Q/i - ^^^3 sin /3) + S,* (Ffcs - ^^3 cos /5)] 



=1 k=\ 



171^171^ 

(C^ - iC[,){K,V*D,^,)^^f{mlml-,mj,^) (32) 



where 



Qij — —7=Ui2Vji 



1 

Si, = -^UaVj2 (33) 



V2 

Fig.(lg) contributes as follows 

Ai^g = L'id sin' ew){UaD,u - K.Ui^D.^i) 

COS t/VV^ ^ ^ o 



^gj^/(m^m^-,m|) (34) 
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where 

1 



4 = -t^at^.i - ^t^;^t^.-2 + 5,j sm' Ow (36) 



The contribution of Fig.(lh) is as follows 



Ai?g = -— E E a, sin Ai^ + a;,Dni) 
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1=1 k=l 



and 



{r]u cos P + r]i* sin P)Y^f{^xl ' ' "^t; ) (3'^) 



^^2] = E E ^fci cos P{atkDai - '^tkDt2i) 



9 



=1 A;=l 



where 



{r]ii cos/? + r;.; sin/3)^^/(m^o, m]j-,m\) (38) 



2 

9mm r-i* I 9m r^* 9 • «n n* ^Q^^ 

—DmDt2j + ^= — AiiAii - ^"^w sm/^AiiAi . (39) 



and 



/ 5'"ib 4 n* n , 9m p, 

\/2mwC0sp vSm^ySin/? 

+^ ^Db2jDt2i + ^= ^^bij Ah - ^"^w cos (3D^^.Dtu (40) 

v2miyC0sp v2m^cosp V2 

Finally the contribution from Fig.(li) is as follows 



= ^ E E T.{Q3siYn + iYis sin /3) + SjsiYi2 + iVis cos P)} 
(f/.iAifc - i^6^s2A2fe)[G'fci(F/2 + ^>/3 cos/3) + G*,(ll2 - cos/3) + iJ,,(Fn + 2^,3 sin/3) 

TTT- — 

+i/4(l^i-zr,3sin/3)]^/(m?^,m^^,m^-)(41) 



3 2 2 
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■(9)_ 9 



V ^ 1=1 s=l k=l 

{KtV:^D,,k)[Gk^{Yl2 + tYis cos/5) + GUY12 - lYi^ cos/5) + Hk^{Yn + lYi^ sin/5) 

Tfl — 

+H:,{Yn - lYi, sin /5)] (< , , mj_ ) (42) 

Summing the contribution from the nine loop diagrams of Fig.(la-li) we find that ARbj 
and ALftjj that appear in Eq. (5) are then given by 



Ait:,,, = E ^< (43) 

n=l 



AL,,, = E (44) 

n=l 

We note that the loop diagram of Fig.(lc) with the interchange Z ^ Hf vanishes in the 
zero external momentum approximation because the vertex is proportional to the external 
momentum, Similarly, the loop diagram of Fig.(lh) with the interchange W~ H~ 
vanishes and the loop diagram of Fig.(li) with the interchange Z ^ vanishes in the 
zero external momentum approximation. 

2.2 Analysis of corrections ARuj and AL^j 

The corrections ARuj and ALuj arise from the nine diagrams of Fig. (2), i.e., the loops 
(a)-(i) of Fig. (2). We label the contribution from the nine diagrams by superscripts 1-9. 
Thus, for example, the contributions of Fig. (2a) are AR[Ij and AL^]^] etc. We now list the 
contributions of the nine loops of Fig. (2). We have 

= 1^ £ KtVj2Dl,,D,,,e^^W,,,mtm^f{ml m^, m?J (45) 



k=l 

A4j = J2iU;iDl^, - K,U*2Dl2k)Dt2iD,2ke-'^'rntrn~J{rnl m?, m? ) (46) 

k=i 

= E E '^KtV^2DlMiDni + A2i) 

k=l 1=1 

iPuDuk + aliD,2k)-^f{ml mjo, m|) (47) 
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where 



and 



k=l 1=1 

rribm^o 

V ^ k=i 1=1 

{Eki{Yi2 + tYi3 cos /?) + E;^{Yi2 - lYis cos /5) + + lYis sin /3) 

- iYis sm/3))(Ci + iC!,)^J{ml m|, m^J 



V^cost/vK ^ o d 



V im^i/ sm p V Zmw sm p 

^^J' = — J^^n ~ ^ ^^^^ dw)DliiDtij + ^ sin^ OwDl^iDt^j] cos (3 
y 2 cos Ow ^ o 6 

\/2mty sin /5 ' 



Cbl = C'b^ cos Xb - C^i sin Xb 
Cbi = (^bl sin Xb + C'ft^ cos 
\f2C^i = Re{hb + + [-Im{hb + 5/^6) sin/3 

+Im{Ah) cos/3]Fz3 + i?e(A/ib)Yi2 
^26-^7 = -/m(/i6 + (5/i6)ya + [-Re{hb + 5hb) sin/3 
+i?e(A/i5) cos P]Yis - Im{Aht)Yi2 



ALg = E E KbU;^Dt,k[Et,,{Yi2 + ^1^3 cos/3) + E*^{Yi2 - lYi^ cos /3) 
+Ffe,(yu + iYi, sin /3) + i^;(yu - iYi, sin /3)] (C,^ - iC^i) {^fiml m| , m^^) 

^< = y E(A/Ai^ + «r,A2.)(i?f, + i?,^)C,cos/3^^/ 



2 3 
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9 1 = 1 ^ 



Ai?lg = 



and 



and 



= 1^ R'i*{atiDni - -ftiDt2i) f{ml rn^o, m^-) 



Ai?;g = - E E [^0- (^1 + ^^fe3 sin (3) + Sij (1^2 + iYks cos 
i=i fe=i 



2 3 

ALjg = E[QM^fei - ^^fc3sin/3) + 5* (^2 - in3Cos/3)] 
«=i fc=i 

rribm- 

{C'k - <u){K,U:,Dm)-^f{ml mj-, mf^J 

cos o 

:^f(mh,m -..my) 



where 



and 



5' ;=i k=l 

Wii sin /3 + r]'i* cos /3) Y^/(m^o , m^- , m? ) 



= EE ^'jkS^^P{(^bkDbii - IhkDm) 
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1=1 k=l 



iv'ii sin P + rfn cos (5) Jg^/ (^x^ > ^i- > ) 
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3 2 2 
V ^ 1=1 s=l k=l 

(KiAifc - KtVs2Dt2k)[Eki{Yi2 + zYi3 cos/3) + - cos/3) + Fk^{Yn + lYis sin/3) 

777. — 

+F:,{Ya - zYis sin /3)] , m^,^ , m^-)(64) 

322 

ALg = E E E{Q;.(>^/i - iYi^ sin/?) + - iYis cos/3)} 

{K,U:^Dak)[Eki{Yi2 + iYi^cos(5) + E;^{Yi2 - iYi^cos(5) + Ffe,(Fn + lYi^ sin/3) 



+i^;(ya - iy,3sin/3)]^/(m|,mi^,mJ-) (65) 
Ai?,,, = E (66) 



AL,,, = E (67) 

n=l 

The loops corresponding to Fig. (2c) and Fig.(2i) where Z <-> vanishes for the same 
reason as discussed earher in Sec. (2.1). Similarly, the loop corresponding to Fig.(2h) with 
W'^ ^ vanishes for the same reason. 



3 The effective Lagrangian for qqi^^j interaction 

We turn now to an analysis of the loop corrections to the squark-quark-neutralino inter- 
action. We begin with the tree level gg^Xj interaction which is given by 

C = gh[KujPR + MujPL\x% 
+gi[KujPR + MujPL]xYti + H.c. (68) 

where 

Mbij = -V2[abjDbu - ^bjP>b2i] 

Mtij = -V2[atjDtu - ltjDt2i\ (69) 
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The loop corrections produce a shift in the couphngs of Eq. (69) as follows 

Ceff = gbiiKbij + AKbij)PR + {Muj + /:^Muj)PL]x% 

+gt[{Kuj + ^Ku,)Pr + {Mu, + ^Mu,)Pl]x% + H.c. (70) 

Thus in this part of the analysis we will calculate the quantities AKbij, AKuj, AMf^j, 
and AMfij from the one loop corrections arising from Figs. (3) and (4) using as in the 
previous anslysis the zero external momentum approximation. 

3.1 Analysis of loop corrections to 66zXj interaction 

Loop corrections to hhiX^j interaction, i.e., AKuj and AM^ij, arise from the nine diagrams 
of Fig. (3). We give now the individual contribution of these nine loops. The contribution 
from Fig. (3a) is 

Ai^g = E e^^W,,,D,,,ial^D;,, - Y,^Dl,,)m,m-J{ml m^, (71) 

-^■"y k=i 

AMg = j- e-'^^D,2iD,2M^D;,, + a,jD;,,)m,m^fiml m|, m^J (72) 

'^^S' k=i 

Fig. (3b) contibutes as follows 
9 1=1 k=i 

"^bm^Ul 2 2 



T^/K'K?'<) (73) 



(2) 2v^ . ^ 



^^hij = E E("M-^Wi ~ lblDb2i){(^blDblk - lblDb2k){PljP>llk + OlbjP>l2k) 

9 1=1 k=i 



-Y^fimlmlo^ml) (74) 
Fig. (3c) makes the following contribution 



= -EE[G'iki(Y,2 + ilia COS /3) +G'*,(ll2 -iYi3 COS /3) 



9 1=1 k=i 



+Hk^{Yn + ^F/3sin/3) + H*,{Yn - lYi^ sin ^C^i + ^C^iM,D,,k + a*D,2k] 



ft 2 2 2 \ /frr-N 

) 
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AMg = - E E \Gh {Yi2 + iYi^ cos (3) + ^^.(yis - iYi^ cos /5) 

1=1 k=i 

+Hki{Yn + ^y^3sin/5) + H*,{Yn - iYissmP)]{C^i - iC^i)[ak,Duk - 7bjDi,2k] 

^& i-/ 2 2 2 \ 
> 



Fig. (3d) contributes as follows 

4 3 

'(4) _ _ 

1=1 k=l 



Axg ^-Y^Y.iQ'ij(^ki + iYk3smP)-Sl^{Yk2 + iYk3COsP)] 



a r> 1711)171^,0 n n n 

= - E E [Qj - ilfes sin /3) - (^2 - iY^z cos 

i=l k=l 



where 



g^ = -=[X3^(X*,.-tan^^X*.)] 



1 



Fig.(3e) conributes as follows 



A ^-(5) 4^251 ,„ 1 . 2a \(Q , *n \'^b'^x'i rf 2 2 2 N 



(5) 472^^ ,,, 1 1 . 2^ V n n 2 2 2x 

" cos^ Ow ^ ^ 2 ~ 3 - 7wA2i)^^^/(mb, m^, m^o) 

/// 1 1 



/// r'"* 



The contribution of Fig.(3f) is 



1=1 

X, r/ 2 2 2\ 
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rritm 



r(6) _ _ 

1=1 

X; /•/ 2 2 2 



167r2 X( 

Fig.(3g) contributes as follows 



/K,m^-,m^-) (84) 



Axg = (85) 



= ^ E ^ii^t^; AH^^/(m^ mj- , m^) (86) 
The contribution from Fig.(3h) is 

2 2 

= - E E ^iA: Sin PiVk.Dtn - KtVk2Dt2i){Vii COS + ri'* sin /?) 



fc=i 1=1 

m 



f{ml^,ml,ml-) (87) 



167r2 



2 2 

= EE^ifc cos/3(ii',t/,*2AiO(% cos/3 + 7^;; sin/?) 

k=l 1=1 



J(m%,mf ,mir-) (88) 
Finally Fig.(3i) gives 

4 2 3 

= E E E[G'fc^(yi2 + iYi3 COS /3) + ^^.(y^s - iYis cos /3) + Hki{Ya + iYi^ sin /3) 

s=l k=l 1=1 

+H:,{Ya - tYi3 sin + alD,2k)[Q'sj{yii + ^>^3sin/5) - S'^^iYi^ + ^^scos/?)] 

'^XS f/_2 _2 _2 



/(m o,mr mj^ )(89) 



= E E i2[Gki{Yi2 + tYis COS /?) + G*,{Yi2 - lY^^ cos (5) + i^fe,(ra + zF^s sin (5) 

s=l k=l 1=1 

+H*^iYn - iYis sin P)]{aMik - lbsDb2k)[Q'*s{Yn - iY^asin/?) - S'*{Yi2 - iY^s cos /3)] 

"^X° //_2 _2 _2 



/(m^o,ml,m^,)(90) 
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The sum of the contribution of the nine diagrams of Fig. (3) gives A/T^jj and AMftjj 
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in) 
bij 
n=l 
9 



AM,,, = Y: AM<J (91) 

n=l 

We note that the diagram corresponding to Fig. (3c) and Fig.(3i) with Z-^Hi, and the 
diagram corresponding to Fig.(3h) with W~ H~ vanish in the zero external momentum 
approximation. 

3.2 Loop corrections to tiiX^ interaction 

Loop corrections to ttiX^ interaction, i.e., AKuj and AMuj, arise from the nine loops of 
Fig. (4). We now give the explicit computation of each loop. Fig. (4a) gives 

^^9 t^i 



= -4z? E e^«^ Ah Ai.«, - lhD;,,)m,m-J{ml m?, m|) (92) 



AMg) = g e-«^ A2. A2.(A*- A*u + c^tjD:,,)rmm-J{ml m?, m|) (93) 

Fig. (4b) gives 

Axg> = j^jZiPtiDm + alMWuDtik + a;,A2.)(«yA*ife - I^D;,,) 

9 1=1 k=i 

^/(m^mJo,m|) (94) 

AMg) = E(««Aii - ltiDt2i){auDnk - luDt2kmjD;,, + atjD;,,) 

9 1=1 k=i 

-^f{mlmlo,ml) (95) 

Fig. (4c) makes the following contribution 

Axg = - E E[^fe^(>i2 + cos/3) + £;4(ll2 - ilia cos/3) 

fi* «=1 A:=l 

+F,,(Fu + tYissmi3) + i^* (Fa - tYissmP)]{Ci + iC^i)[(3t,Dnk + «t,A2fc] 



^f{mlm]j^,ml) (96) 
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3 2 



AMg) = -Y.Y.[Eki(Yi2 + iYisCOsP) + E*,(Yi2-iYisCOsP) 

9 1=1 k=i 

+Fu(Yn + iYasmP) + FHYu - tVasynmcS - iC^)[at,Dat - 7yA2»] 

^/(m?,m|,„m|) (97) 

Fig.(4d) gives 

Axg = -EE[Q;,(nl+^n3sm/3) -5;^.(n2+in3cos/3)] 



(Cf, + iCf,){PuDni + alDni)^^f{mt, m^o, m^J (98) 



AMg) = - E E [QiKni - ^n3 sin /3) - 5^.^ (^2 - lY^z cos 



4 3 

EI 

1=1 fc=i 



(Qi - tCf,){auDai - ltiDt2i)-^ f{mlmlo,m\) (99) 



Fig.(4e) gives 



= ;;^E^;';(^in'^^)(A/AH + a:,A20^^/(m^m|^ (100) 



"^^^ ^ ^^^^^^'^^ - ^sin^MKAn -7«A20^^/(m?,m|,m^o) (101) 
Fig.(4f) gives 

= E(^w - ^wKV^aAi^ - K,VnDt2i%isiTi(3 
1=1 

rribm + 

-^f{mlml^,ml^) (102) 



AMg) = -E(i^6t + 5w*)(^;>t^r2Aii)e.7cos/3 
z=i 



Fig.(4g) gives 



'^''^''^ f{mlml+,ml+) (103) 



Axg = (104) 
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AMg) = ^ E i?;. W2Ai.^^^/(m^, mj+, m^) (105) 
Fig.(4h) makes the following contribution 

2 2 

Axg = - E E ^Ifc cos (5{UkiD,u - K,Uu2Dm) {v'u sin P + Vu cos /?) 

A;=l Z=l 



2 2 

AMg.^ = E E & /3(i^tV^;2^6i/) (^H sin P + r;* cos (3) 

k=l 1=1 

^2 ^2 



Finally Fig.(4i) gives 

AkS] = E E Yl[Eki{Yi2 + iYis cos /3) + E*,{Yi2 - iYi^ cos /3) + Ffc,(Y,i + zY^g sin /3) 

s=l A:=l ;=1 

+F*,{Yn - lYi, sin P)]{Pt,Dnk + <A2fc)[g;,(>ii + lYi^ sin f3) - S^^iYi^ + lYi, cos/3)] 

^/Ko'<' 0(108) 

AMg) = E E j2[Eki{Yi2 + ^l/3 cos /3) + E*,{Yi2 - lYi^ cos (5) + + zF,3 sin (5) 

s=i it=i ;=i 

+i^;(yu - il/3sin/3)](at,Aiik - 7t.A2fe)[Q;:(l/i - ills sin /3) - 5;:(Y;2 - iF/aCOS/?)] 



267^2-' ^ X^' til 

The sum of contributions above give AKbij and AMinj so that 



n=l 

am,,,=x:a^s^ (110) 



n=l 



As in the previous analysis, the contribution from Fig. (4c) with the interchange Z <-> 
vanishes in the zero external momentum approximation since the vertex is proportional to 
the external momentum. Similarly, the contribution from Fig.(4h) with the interchange 
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^ H vanishes and Fig.(4i) with the interchange Z ^ Hi vanishes in the zero 
external momentum approximation for the same reason. We also note that loops where 
one of the internal lines is a gluon line also vanishes in the zero external momentum 
approximation since the squark-gluon interaction gives a vertex of —igsip + p')'^ which is 
of course dependent on the momenta. 

4 Loop corrected squark decays into charginos and 
neutralinos 

Eqs.(5) and (70) give the loop corrected effective Lagrangian for qq^xf ^-nd qqiX^ inter- 
actions. Next we use this loop corrected Lagrangian to compute the decay widths of the 
third generation squarks into charginos and neutralinos. Specifically we will analyze the 
following decays 

bi^t + xj 
ii^b + xt 
bi^b + x] 

ii^t + x'j (111) 

To make the analysis more compact we begin by writing both Eqs.(5) and (70) in the 
following form 

C = /(4 + Bf^^s)fjqi + H.c. (112) 

where / takes on the values (t, b) and fj stands for xf, X° while q^ can be bi, ii. The decay 
width r{qi fjf) is given by 

{^(l^r + m){ml - m] - m)) - l(|i^gf - \BZ?)2m,m,} (113) 

The co-efficients Bf^ and B^^ contain the loop corrections and depend on the CP phases.Thus, 
for example, the process &i — > + ^ gives the co-efficients 
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= 2^^^j ~ ^bij + ^Rbij — ^Lhij] (114) 
where Rbij, Luj, ALuj and ARuj are defined by Eqs.(2), (43) (44). 

5 Numerical Analysis and Size of Effects 

In this section we discuss in a quantitative fashion the size of loop effects on the decay 
widths of the squarks into chargino and neutralinos. The analysis of Sees. (2-4) is quite 
general and valid for the minimal supersymmetric standard model (MSSM). For the sake 
of numerical analysis we will limit the parameter space by working within the framework 
of the SUGRA models[23]. Specifically within the framework of the extended mSUGRA 
model including CP phases, we take as our parameter space at the grand unification 
scale to be the following: the universal scalar mass mo, the universal gaugino mass mi/2, 
the universal trilinear coupling \Aq\, the ratio of the Higgs vacuum expectation values 
tan/3 —< H2 > / < Hi > where H2 gives mass to the up quarks and Hi gives mass 
to the down quarks and the leptons. In addition, we take for CP phases the following: 
the phase of the Higgs mixing parameter /i so that /i — |/i|e*^'', the phase a^o of the 
trilinear couphng where = | Ao|e'"^o , and the phases (i=l,2,3) of the SU (3)c, SU {2)l 
and U{1)y gauginos, so that rhi — |mi|e'^* (i=l,2,3) where rui (i=l,2,3) are the SU{3)c, 
SU{2)l and U{1)y gaugino masses. We note that not all the phases are independent and 
only certain combinations of them appear in the analysis [5]. In the numerical analysis we 
compute the loop corrections and also analyze their dependence on the phases. 

In Fig. 5 (a) we give a plot of the decay width of the heavy stop (ti) into light and heavy 
chargino, xt and xt } i-^-i ^ P^ot of r(ti &Xi^2) ^-s a function of aA„- The plots arc given 
with the analysis done at the tree level and at the level of the cfTcctivc Lagrangian including 
loop corrections. The analysis shows that the loop effects can produce a correction of as 
much as 25% to the tree level values. Further, the analysis of Fig.5(a) shows that the 
dependence on aA„ is quite significant and both the tree and the loop corrections are 
affected by it. From Fig. 5(a) one finds that the variation with in the range (0, tt) can 
be as much as 40-50%. In Fig. 5(b) a similar plot is given for the decay width T{ii txl 4) 
as a function of a^o- Here one finds that the loop corrections can be as much as 20% and 
further that the variations with ct^g can be as much as 25-30%. The effect of on the 
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decay width arises from two sources: (1) enters the off diagonal elements of the squark 
mass^ matrix. So it affects the squark masses that enter in the decay width. In fact, the 
modification of the squark masses due to a^o can be large enough that a decay channel 
may close or open as ce^g is varied. This phenomenon will be illustrated explicitly later. 
This type of effect appears both at the tree and at the loop level. (2) The matrix Dq that 
diagonalizes the squark matrix is sensitive to variations of and this variation again 
affects both the tree and the loop level analysis. Thus at the tree level the couplings Rqij, 
Lgij, Kqij and Mq^j depend on and similarly at the loop level the couplings ARqij, 
ALqij, A.Kgij and AMqij also depend on a^g. An important phenomena related to the 
dependence on is that the effects are strongly dependent on the quark mass. This 
is so because phases enter in the squark mass^ matrix via the off diagonal terms in a 
prominent way and these off diagonal terms are proportional to the quark mass. Because 
of this, the sensitivity of the stop decay widths to ckaq is far greater than the sensitivity 
of the sbottom decay width. The loop corrections are bigger in the case of the stop decay 
than for the sbottom case due to the relative difference of their Yukawa couphngs. For 
this reason in our numerical analysis we will focus mostly on the effects of phases on stop 
decays. 

Fig. 6 (a) is a repeat of Fig. 5 (a) with a plot of the light stop (^2) decay width into 
charginos, i.e., r(t2 — ^ a function of a^^. Here one finds that while the loop 

corrections are comparable to the case of Fig. 5(a), the variations of the decay width is 
more strongly dependent on «a„ in this part of the parameter space. Fig.6(b) gives an 
analysis similar to that of Fig. 5(b) where plots are given for the decay width T{t2 tx^ 4) 
as a function of a^g. Here one finds that the loop corrections can be as much as 30%. 
Further, one finds that the variations with a^^, are now much stronger than in the case 
of Fig. 5(b). Thus the effect of aA„ is large enough that for values of aAo > 1-3 (radian) 
the decays into Xs^ Xa ^t^^ closed. The reason for this is purely kinematical, in that the 
mass of ^2 is strongly dependent on oaq and varies strongly with aAo ^^^^^ below the 
kinematical limit to allow for the decay into X3, xt for values of aAo — l-^- In Fig-7(a) a 
plot is given of the decay width r{ti —>■ bx'^, tx°) (where we summed over the final states 
of charginos and neutralinos) both at the tree level and at the loop level as a function of 
aAo- 

The analysis of Fig.7(b) is similar to that of Fig.7(a) except that one is looking at 
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the decay width of t2- The discontinuities in Fig. 7(b) are kinematical and arise from the 
closing of some of the neutrahno final states. The analysis of Fig. 8 (a) is similar to that 
of Fig. 7(a) while the analysis of Fig. 8(b) is similar to that of Fig. 7(b) except that the 
plots are made as a function of 6^. It is interesting to observe that the dependence of 
the stop widths on 6*^ in Fig. 8(b) appears to be relatively weaker. This arises because we 
are summing over the chargino and neutralino final states. Thus, for example, the decay 
width r(ti — >■ bxi) increases with for the parameters of Fig.8(a) while r(ii — > bxt) 
decreases. This results in the sum r{ti — > bxi, bxt) having only a weak dependence on 
9^. An analysis similar to that of Figs.7(a)-7(b) but as a function of ^3 is carried out 
in Figs.9(a)-9(b). One important new feature of the decay widths here is that the ^3 
dependence of the widths at the tree level is absent while the loop corrected widths show 
a dependence on ^3. Here one finds that the loop corrections are typically of size 10-15% 
while the overall variation with ^3 can be as large as 20%. Typically, the loop correction 
to the sbottom decays are small and the dependence on phases is also relatively small. 
This is exhibited in Fig. 10 where the decay width r(bi — > txi) is plotted. Here one finds 
that the loop effects are essentially negligible while the variations of the decay width with 
ccao is also essentially negligible. The reasons for this weak dependence on the phase and 
the smallness of loop corrections have already been explained on an analytical basis at 
the end of the second paragraph of this section. Here we see that the reasoning presented 
there is borne out by the numerical analysis. Thus the largest loop corrections as well as 
the largest variations with phases arise only for the decay of the stops. 

The experimental upper limits of the electric dipole moments are[8, 9, 10]: \de\ < 4.3 x 
10~^''ecm, \dn\ < 6.5 x 10~^^ecm and \dHg\ < 9.0 x 10~^^ecm. The last constraint for Hg^^^ 
could be transformed into a constraint on a specific combination of the chromoelectric 
dipole moments of u, d and s quarks[6], Cug = ~ du ~ 0.012(if | < 3.0 x lO^^^cm. 
These constraints are satisfied by the cancellation mechanism in the numerical analysis 
presented above as follows: In figure 5, 6 and 7 the constraints arc satisfied for the inputs 
tan/? = 40, mo = 300 GeV, mi/2 = 300 GeV, ^1 = 0.5 (radian), ^2 = -66 (radian), 
^3 = .63 (radian), = 2.5 (radian), a^o — 1-0 (radian) and |Ao| = 1. At this point 
we have |4| = 1-88 x lO^^^ecm, \dn\ = 1.79 x lO'^^ecm and Cng = 8.99 x lO'^W. In 
figures 8, 9 and 10 they are satisfied for the inputs tan P — 45, mo = 400 GeV, mi/2 = 400 
GeV, ^1 = 0.6 (radian), ^2 = -65 (radian), ^3 = .65 (radian), 9^ — 2.5 (radian), — 2.0 
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(radian) and \Ao\ = 1. At this point we have \de\ = 3.94x10 ^''ecm, \dn\ = 9.21x10 ^''ecm 
and Cffg = 3.86 x lO-^^cm. 

6 Conclusion 

In this paper we have analyzed supersymmetric one loop corrections to the squark-quark- 
chargino and to the squark-quark-neutralino couplings. The analysis involves the exchange 
of the gluino, chargino, neutrahno, W, Z, charged Higgs and neutral Higgs. With the 
above analysis the one loop effective Lagrangian for these interactions was derived. The 
full CP dependence arising from the soft CP parameters was taken into account in the 
analysis. The effective Lagrangian was then used to obtain the decay of the squarks into 
charginos and neutralinos at the one loop order. A detailed numerical analysis within 
extended SUGRA model was then carried out to study the size of the loop effects and 
also to study the effect of CP phases on the decay widths of the squarks into charginos 
and neutrahnos. The analysis exhibits that the loop corrections to the decays widths 
of the stops can be very substantial, i.e., as much as 30% or more. Further, the phase 
dependence of the decay width is found to be very strong producing a variation of as much 
as 40-50% or more. The phases enter in the decay widths in two ways; in modifying the 
stop-bottom-chargino, and the stop-top-neutralino couplings and in modifying the stop, 
chargino and neutralino masses. In some cases the effect of phases is large enough to 
open or close a decay channel. However, a similar analysis for the decay of the sbottoms 
shows the effect of loops as well as the effect of CP phases to be much smaller. The one 
loop effective Lagrangian derived in this paper would be useful in the analysis of squark 
decays at colliders and in connecting experimental data with the underlying theoretical 
schemes such as supergravity and string based models. 
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Appendix A 

For completeness we give below the loop corrections to the Yukawa couplings A/i^, Shj, etc 
that appear in Sec. 2. A derivation of these results can be found in Ref.[13] and Ref.[15]. 
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2 2 2a 
i=ij=i '^^ 

2 2 2 fj2 + 

- E E E 9'e:^{V,\D;,, - K,V,*,D;,,}{K,U;,Dnj)j^J{ml,,rnl m? ) 

i=l j=l k=l 

2 2 2 772 +m + 

- E E E 9"C^Ay^lD*a, - i^*V;;i^;2j(^.^2Ai.)^^/(m|, m^,, 

i=l j=l fc=l 



2 2 4 



+ E E E 2Q,{«6fcAi, - 7bkD,2j}{f3;,Dl,, + a,,Dt,,}-^J{mlo ml,ml) 
i=i j=i k=i 



4 4 2 



+ E E E 2r.,{a,,- Aife - 76.- A2a{/96* i^Hfe + aMDl,,}^^f{ml, mjo, mjo) (115) 
i=i j=i k=i ' ^ 



v/2 sin/3^2Mw^ 



% -<5^jCos/?-i^*,• 



yJ 



(116) 



r 



^7 r "^x? 



V2 2sin/?^2Miy 



%-g,;cos/?-i?,;] 



(117) 



~ 2M, 



9(^ij - T^ixUgx;, - g'xi^) + (i ^ j)] 



R 



2Mi 



(118) 



2 2 

i=ij=i -^^ 

2 2 2 172 , 

-EE E - i^*V^;2A*2j(^^.f/:2Ai,)Y^/(mJ,, m|, m|) 

2 2 4 ^0 

+ E E E 2i^,.{«6fcAi, - 76./^62,}{/34^M. + (^bkDt2^}j^f{mlo,ml,ml) (119) 
1=1 1=1 fe=i 



i=ij=i 
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2 2 2 fJl : 

- E E E 9'h:.{U;,DI,, - K,U;,Dl,,}(K,V:,D,,,)^J(rnl,, ml m|) 

i=l j=l k=l 

2 2 2 m+m + 

E E E 9'K*,{U:,Dl,, - K,U:,Dl,,}iK,V;M^^firnl, mj,, m^,) 

i=l j=l k=l 

2 2 4 Q 

+ E E E 2i^H«*^^*y - mA2,-}{A;£'a, + ^tkD:,,}j^J{mlo,ml,ml) 



1 j=i fe=i 

4 4 2 m^om^o 



+ E E E 2A^{a,,-Ai. - 7*, A2j{A*.A*u + rnjo, mjo) (120) 

i=i j=i fe=i ' 

where 

K,^ = -v^^Q.i (121) 

= -^g;; (122) 

2 2 2a 

= - E E -^e-^^'rn-,E,,D;,,Dt,,f{ml m|, m|) 

2 2 2 ^ , 

- E E E 9'G,dU:,D;,^ - K,U*M{K,V,%D,,,)^J{ml,,ml^^^^ 

i=l j=l k=l 

2 2 4 „^ „ 

+ E E E 2^.^{«tfcAi, - -ftkDt2jm,Dl,^ + at,D;,^}:^f{mlo,ml,ml) (123) 
1=1 j=i fe=i 



-A/i, = - E E ^e-^^W,,jD:,,rt;,m-J{ml m|, m? ) 
i=ii=i 

2 2 4 

+ E E E ^ii^^bkDbij - lbkDb2jWtkDni + <^tkDl^d 

1=1 j=i k=i 

^vO 2 2 4 m^om- 

^/(mjo , m|, m|.) + ^ E E V2gCki^^:^[-K,U:,Dnj{P:,D;,^ + atkO;,^) 
f{ml.,ml-,mlo) + {a.kDnj - i,kD,2j)iU;,Dl^^ - K,U:,Dl,^)f{ml,ml-,mlo)] (124) 



-Sh = E E ^e-^«3^^^^.^*^^^^.^^,y(^2^ 

2 2 4 ^ Q 

- E E E 2^M«6/c/^6i, - %kD,2jMkD;u + «t/iA*2^)T^/(rnx2' ^' <) (125) 
i=i j=i k=i 
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-Ah, = - E E ^e-^''D;,,D,,jrj'*rn,f{rnl m? , m|) 

+ E E E '^nCMtkDtij - -ftkDt2j){PtkDlu + abkDl2i)-7^f{mlo,ml,ml) 
i=ij=ik=i -"-"^ 

2 2 4 fYj^ Q <yj^ _ 

+ E E E ^9Ck.^^[-KtV^D,^j{Pt,D;,, + a,,Dl,,) 

i=i j=i k=i -"-"^ 

f{ml,ml-,mlo) + (a^fcAi, - 7tfeA2,)(^aAy - ^^^i^A^O/K- ^o)] (126) 

^ ^ 2q; 

-^t = E E -V^(^~'^'D;^i^t2jVijmgf{ml, m? , m|.) 

2 2 4 y 

- E E E 2%(at;t Ai,- - ltkDt2jmkDlu + «6feA*2j ^) (127) 

i=l ;?=1 fc=l 
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Figure 1: List of one loop graphs that contribute to the bitXj coupUngs arising from the 
exchange of the gluino, charginos, neutrahnos, W, Z, charged Higgs and neutral Higgs. 




Figure 2: List of one loop graphs that contribute to the tibxj couplings arising from the 
exchange of the gluino, charginos, neutralinos, W, Z, charged Higgs and neutral Higgs 
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Figure 3: List of one loop graphs that contribute to the b^bx^ coupUngs arising from the 
exchange of the gluino, charginos, neutrahnos, W, Z, charged Higgs and neutral Higgs 
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X? * X? * x5 
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Figure 4: List of one loop graphs that contribute to the UiXj couphngs arising from the 
exchange of the gluino, charginos, neutralinos, W, Z, charged Higgs and neutral Higgs. 
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O 0.5 1 1.5 2 2.5 3 

(a) Plot of the decay width T{ii &Xi^2) ^ ^ function of uaq- The sohd hnes correspond to 
analysis at the tree level while the long-dashed lines include loop corrections. The inputs for the 
thin lines is tan/3 = 40, mo = 300 GeV, mi/2 = 300 GeV, = 0.5 (radian), ^2 = -66 (radian), 
^3 = .63 (radian), = 2.5 (radian), and \Ao\ = 1. The thick lines are for xt decay and the 
thin lines are for xt decay. 
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O 0.5 1 1.5 2 2.5 3 

(b) Plot of the decay width r{ii txs 4) as a function of a^o- The solid lines correspond to 
analysis at the tree level while the long-dashed lines include loop corrections. The inputs for the 
thin lines is tan/3 = 40, mo = 300 GeV, mi/2 = 300 GeV, ^1 = 0.5 (radian), ^2 = -66 (radian), 
^3 = .63 (radian), 9^ = 2.5 (radian), and \Ao\ = l.The thick lines are for x^ decay and the thin 
lines are for X3 decay. 
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(a) Plot of the decay width T{i2 &xf,2) ^ ^ function of uaq- The sohd hnes correspond to 
analysis at the tree level while the long-dashed lines include loop corrections. The inputs for the 
thin lines is tan/3 = 40, mo = 300 GeV, mi/2 = 300 GeV, = 0.5 (radian), ^2 = -66 (radian), 
^3 = .63 (radian), = 2.5 (radian), and \Ao\ = 1. The thick lines are for xt decay and the 
thin lines are for xt decay. 




(b) Plot of the decay width r(t2 txs 4) as a function of a^o- The solid lines correspond to 
analysis at the tree level while the long-dashed lines include loop corrections. The inputs for the 
thin lines is tan/3 = 40, mo = 300 GeV, mi/2 = 300 GeV, ^1 = 0.5 (radian), ^2 = -66 (radian), 
^3 = .63 (radian), 6^ = 2.5 (radian), and \Ao\ = 1. The thick lines are for Xa decay and the 
thin lines are for Xs decay. 
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(a) Plot of the decay width r{ti bx^,tx'^) as a function of aAo- The sohd hnes correspond 
to analysis at the tree level while the long-dashed lines include loop corrections. The inputs 
for the thin lines is tan/3 = 40, mo = 300 GeV, mi/2 = 300 GeV, = 0.5 (radian), ,^2 = -66 
(radian), ^3 = .63 (radian), 6^ = 2.5 (radian), and \Aq\ = 1. The thick lines are for the sum 
over the neutralino final states and the thin lines are for the sum over the chargino final states. 
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(b) Plot of the decay width T{i2 bx^,tx^) as a function of a^„. The solid lines correspond 
to analysis at the tree level while the long-dashed lines include loop corrections. The inputs 
for the thin lines is tan/3 = 40, mo = 300 GeV, mi/2 = 300 GeV, ^1 = 0.5 (radian), ^2 = .66 
(radian), ^3 = .63 (radian), 6^ = 2.5 (radian), and \Ao\ = 1. The thick lines are for the sum 
over the neutralino final states and the thin lines are for the sum over the chargino final states. 
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(a) Plot of the decay width T{ii bx'^,tx^) as a function of 6^. The sohd lines correspond 
to analysis at the tree level while the long-dashed lines include loop corrections. The inputs 
for the thin lines is tan/3 = 45, mo = 400 GeV, mi/2 = 400 GeV, = 0.6 (radian), ^2 = -65 
(radian), ^3 = .65 (radian), a^o = 2 (radian), and \Aq\ = 1. The thick lines are for the sum 
over the neutralino final states and the thin lines are for the sum over the chargino final states. 




(b) Plot of the decay width r(<2 &X^jiX°) as a function of 9^. The solid lines correspond 
to analysis at the tree level while the long-dashed lines include loop corrections. The inputs 
for the thin lines is tan/3 = 45, mo = 400 GeV, mi/2 = 400 GeV, ^1 = 0.6 (radian), ^2 = .65 
(radian), ^3 = .65 (radian), = 2 (radian), and |Ao| = 1. The thick lines are for the sum 
over the neutralino final states and the thin lines are for the sum over the chargino final states. 



Figure 8: 
34 




(a) Plot of the decay width T{ii bx'^,tx^) as a function of ^3. The sohd hnes correspond 
to analysis at the tree level while the long-dashed lines include loop corrections. The inputs 
for the thin lines is tan/3 = 45, mo = 400 GeV, mi/2 = 400 GeV, ^1 = 0.6 (radian), ^2 = -65 
(radian), 6^ = 2.5 (radian), a^o = 2 (radian), and \Ao\ = 1. The thick lines are for the sum 
over the neutralino final states and the thin lines are for the sum over the chargino final states. 
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(b) Plot of the decay width T{i2 bx'^,tx°) as a function of ^3. The solid hnes correspond 
to analysis at the tree level while the long-dashed lines include loop corrections. The inputs 
for the thin lines is tan/3 = 45, mo = 400 GeV, mi/2 = 400 GeV, ^1 = 0.6 (radian), ^2 = .65 
(radian), 6^ = 2.5 (radian), uaq = 2 (radian), and \Ao\ = 1. The thick lines are for the sum 
over the neutralino final states and the thin lines are for the sum over the chargino final states. 
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Figure 10: Plot of the decay width T{bi — » txi) as a function of The sohd hues 

correspond to analysis at the tree level while the long-dashed lines include loop corrections. 
The inputs for the thin lines is tan/? = 45, mo = 400 GeV, mi/2 = 400 GeV, = 0.6 
(radian), ,^2 = -65 (radian), ,^3 = .65 (radian), = 2.5 (radian), and \Ao\ = 1. The input 
for the thick lines is the same as for the thin lines except that ^2 = -5 (radian). 
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